J Wood Sci (2011) 57:203-207 
DOI 10.1007/sl0086-010-l 157-7 


© The Japan Wood Research Society 2011 


ORIGINAL ARTICLE 


Noboru Nakamura • Kazuhiko Fujita 

Simulation method to generate the strength of glulam using correlated 
random variables 


Received: May 24, 2010 / Accepted: September 29, 2010 / Published online: March 17, 2011 


Abstract Many reports have been published about design¬ 
ing the strength of glulam using simulation methods. In 
simulation methods, one of the most important problems is 
how to deal with correlations among strength factors, i.e., 
modulus of elasticity (MOE), modulus of rupture (MOR), 
tensile strength (ay), and compression strength (a c ). For 
example, in the case that the MOR criteria of glulam is aj 
fni + a hJfbi -1 (where a„, and a bi are the axial stress and the 
bending stress of the i-th lamina respectively, and f ni and f bi 
are the axial strength and the bending strength of the Z-th 
lamina respectively), a correlation betweenand f bi exists. 
How can we account for this correlation when calculating 
the strength of glulam, bearing in mind that it is very diffi¬ 
cult to measure the correlation coefficients among MOR, 
a T , and a c ? We developed a method by which these prob¬ 
lems could be solved, and, using random variables gener¬ 
ated by this method, the strengths of glulam were simulated. 
The simulated values were almost the same as the experi¬ 
mental values. The results indicated the usefulness of the 
method. 

Key words Simulation method • Glulam strength • Correla¬ 
tion coefficients ■ Lamina strength 


Introduction 

Several simulation methods 1-7 for designing the strength of 
glulam have been proposed. In the simulation method, an 
important problem is how to deal with correlations among 
strength factors such as modulus of elasticity (MOE), 
modulus of rupture (MOR), tensile strength (a T ), and com- 
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pression strength (a c ), because some of them are involved 
in the strength criteria of glulam. For instance, in the case 
of two factors, MOE and a x , the simulation method to gen¬ 
erate them using correlated random variables was shown in 
Taylor and Bender, 8 but no method has been put forward 
to incorporate additional factors. The reason might be that 
it is very difficult to measure the correlation coefficients 
among MOR, a T , and a c . Even if the matched specimen 
method is used, we cannot measure the strength factors of 
the exact same specimens. 

Therefore, we developed a method by which the correla¬ 
tion coefficients among the strength factors could be calcu¬ 
lated. Using this method, the strengths of glulam were 
simulated. To validate the method, experiments were carried 
out to determine MOE, MOR, a T , and a c of glulam; the 
results are given in this report. 


Materials and methods 

Specimens and experiments for laminae 

MOE, MOR, a T , and a c were measured for Sugi (Japanese 
cedar) and Douglas-fir laminae. Before tension and com¬ 
pression testing, MOE values were measured by bending 
tests, as shown in Fig. 1. When a T and a c were measured, 
Young’s modulus in tension (E T ) and compression ( E c ) 
were also measured. E c and E T were measured by displace¬ 
ment transducers fixed on four sides of the specimens. The 
grades of Sugi laminae were L30-L90 and those of Douglas 
fir laminae were L90-L160. The grades for laminae were 
determined according to the minimum MOE measured 
using a continuous MOE measuring machine. Compatibility 
conditions for each lamina grade in the Japanese Agricul¬ 
tural Standard (JAS) for glulam are shown in Table 1. In 
this table, 5% values are calculated under the assumption 
that the strength is normally distributed. 

The experimental setup is shown in Fig. l.The width and 
thickness of laminae were 105 mm and 30 mm, respectively, 
and each lamina had a single finger joint. MOE values of 
each specimen for tensile and compression strengths were 
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measured before strength testing and the grade was decided. 
The number of specimens for each grade was about 30. 

Correlation coefficients between strength factors of laminae 

Each lamina broke at the finger joint. The results of the 
strength experiments are summarized in Table 2; in this 
table, the best-fit probability distributions are also shown. 
Kolmogorov-Smirnov examination was used to fit the 
strength data to the probability distribution. 

The relations between MOE and strengths in Douglas fir 
and Sugi are shown in Figs. 2 and 3, respectively. Using these 
regression analyses, correlation coefficients between MOE 
and the strengths were calculated. The results are shown in 
Tables 3 and 4. 


used. Therefore, we calculated the values by the following 
method. First, a combination of uniform random variables 
from 0 to 1 with no correlation was generated by the Mer- 
senne twister method. 9 The number of generated random 
variables for each strength factor was 1000. Next, these 
variables were transformed into a combination of standard 
normal variables with no correlation. Third, they were trans¬ 
formed into a combination of standard normal variables 
with correlation. Fourth, they were further transformed into 
non-normal variables with correlation by Rosenblatt trans¬ 
formation. Using this method, correlation coefficients 
among the strength factors such as MOR, cj t , and a c were 
calculated. The results are shown in Tables 5 and 6. We can 
see that the values measured by experiments are almost 
same as those generated by this method. 


Calculation of correlation coefficients among strength 
factors other than those measured by experiment 

In Tables 3 and 4, correlation coefficients could not be mea¬ 
sured rigorously, even though matched specimens were 
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Validation of the method 

In order to validate the method, using all correlation coef¬ 
ficients in Tables 5 and 6, 1000 strength values for each 
factor were generated and correlation coefficients were cal¬ 
culated by the above-mentioned method. The results are 
shown in Tables 7 and 8. Again, we can see that the values 
in Table 7 and 8 are almost the same as those in Table 5 and 
6. These results indicate the usefulness of the method 
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Fig. 1. Setups for various experiments on lamina 


Table 1. Compatibility conditions for modulus of elasticity (MOE) in 


Japanese Agricultural Standard 


Grade 

Minimum MOE 
(kN/mm 2 ) 

Bending strength 
(N/mm 2 ) 

Tensile strength 
(N/mm 2 ) 

Mean 

5% value 

Mean 

5% value 

L200 

20.0 

81.0 

61.0 

48.0 

36.0 

L180 

18.0 

72.0 

54.0 

42.5 

32.0 

L160 

16.0 

63.0 

47.5 

37.5 

28.0 

L140 

14.0 

54.0 

40.5 

32.0 

24.0 

L125 

12.5 

48.5 

36.5 

28.5 

21.5 

L110 

11.0 

45.0 

34.0 

26.5 

20.0 

L100 

10.0 

42.0 

31.5 

24.5 

18.5 

L90 

9.0 

39.0 

29.5 

23.5 

17.5 

L80 

8.0 

36.0 

27.0 

21.5 

16.0 

L70 

7.0 

33.0 

25.0 

20.0 

15.0 

L60 

6.0 

30.0 

22.5 

18.0 

13.5 

L50 

5.0 

27.0 

20.5 

16.5 

12.0 

L40 

4.0 

24.0 

18.0 

14.5 

10.5 

L30 

3.0 

21.0 

16.0 

12.5 

9.5 


Table 2. Strength properties of laminae 


Species 

Strength factor 

MOE 

(kN/mm 2 ) 

MOR 

(N/mm 2 ) 

Oj 

(N/mm 2 ) 

etc 

(N/mm 2 ) 

Douglas fir 

No. of specimens 

420 

150 

145 

125 


Mean 

13.68 

58.61 

40.35 

47.90 


Coefficient of variation 

0.17 

0.21 

0.21 

0.13 


Probability distribution function 

3P-Weibull 

Normal 

Normal 

Log-normal 

Sugi 

No. of specimens 

498 

170 

158 

170 


Mean 

5.49 

34.57 

23.76 

31.73 


Coefficient of variation 

0.34 

0.25 

0.28 

0.17 


Probability distribution function 

Normal 

Log-normal 

Log-normal 

Log-normal 


MOR, modulus of rupture; cr x , tensile strength; <7 C , compression strength; 3P, three-parameter 
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Fig. 2. Relations between 
modulus of elasticity ( MOE ) and 
modulus of rupture ( MOR), 
tensile strength ( a T ), and 
compression strength (<7 C ) in 
Douglas fir 
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Fig. 3. Relations between MOE 
and various strength indicators 
in sugi 



Table 3. Calculated correlation coefficients for Douglas-fir laminae 



MOE 

MOR 

E t 

(Jf 

E c 

Ob 

MOE 

1 

0.741 

0.881 

0.585 

0.810 

0.893 

MOR 


1 

? 

? 

? 

? 

E t 



1 

0.503 

? 

? 

(J-y 




1 

? 

? 

E c 


Sym. 


1 

0.745 

Ob 






1 


E t , Young's modulus in tension; E c , Young's modulus in compression; 
?, unknown 


Table 4. Calculated correlation coefficients for sugi laminae 



MOE 

MOR 

E t 

Or 

E c 

Oc 

MOE 

1 

0.628 

0.977 

0.731 

0.943 

0.850 

MOR 


1 

? 

? 

? 

? 

E t 



i 

0.718 

9 

? 

Ox 




1 

? 

? 

E c 


Sym. 


i 

0.805 

o c 






1 


Table 6. Correlation coefficients for sugi laminae by the proposed 
method 



MOE 

MOR 

E t 

0"x 

E c 

o c 

MOE 

1 

0.638 

0.979 

0.754 

0.951 

0.855 

MOR 


1 

0.619 

0.476 

0.596 

0.525 

E t 



1 

0.732 

0.929 

0.838 

0"x 




1 

0.718 

0.652 

£ c 


Sym. 


1 

0.810 

o c 






1 


Table 7. Correlation coefficients for Douglas-fir laminae using corre¬ 
lated random variables 

MOE 

MOR 

E t 

< 7 t 

E c 

o c 

MOE 1 

0.744 

0.878 

0.566 

0.807 

0.880 

MOR 

1 

0.673 

0.434 

0.617 

0.662 

E t 


1 

0.488 

0.714 

0.786 

< 7 t 



1 

0.454 

0.502 

E c 

Sym. 


1 

0.742 

Ob 





1 


Table 5. Correlation coefficients for Douglas-fir laminae by the pro¬ 
posed method 



MOE 

MOR 

E t 

Oj 

E c 

o c 

MOE 

1 

0.749 

0.882 

0.583 

0.809 

0.892 

MOR 


1 

0.669 

0.464 

0.600 

0.674 

E t 



1 

0.508 

0.731 

0.797 

<7 t 




1 

0.459 

0.517 

E c 


Sym. 


1 

0.736 

O c 






1 


Table 8. Correlation coefficients for sugi laminae using correlated 
random variables 



MOE 

MOR 

E t 

<J T 

E c 

o c 

MOE 

1 

0.647 

0.976 

0.760 

0.946 

0.853 

MOR 


1 

0.638 

0.480 

0.605 

0.524 

E t 



1 

0.742 

0.932 

0.829 

0"x 




1 

0.735 

0.661 

E c 


Sym. 


1 

0.805 

o c 






1 
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Fig. 4. Specimens for glulam experiments. The width of all glulam 
specimens was 105 mm. The thickness of type A and C specimens was 
150 mm and the thickness of type B and D specimens was 300 mm. 
Shaded laminae, Douglas-fir; white laminae, sugi 
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Fig. 5. Setups for the glulam experiments 


proposed here. In this method, all correlation coefficients 
were transformed according to the best-fit probability 
distributions. 10 

Simulation of glulam strengths using correlated 
random variables 

Using random variables generated by the above-mentioned 
method, glulam strengths were generated and compared to 
those of experiments. 

Specimens and experiments of glulam 

Glulam specimens are shown in Fig. 4, in which the combi¬ 
nations of laminae are also shown. Bending, tensile, and 
compression strengths were measured in about six speci¬ 
mens for each experiment. Bending experiments were con¬ 
ducted according to JAS, but compression experiments 
were conducted according to Testing and Evaluation 
Method for Full-size Structural Timber Strength, 11 as shown 
in Fig. 5. In contrast, tension experiments were not con¬ 
ducted to the specified standard, because the length of 
specimens was not as long as shown in Fig. 5. 


Results and discussion 

An example of the relationships between strength and 
cumulative probability for type A and B are shown in Fig. 
6, and also type C and D in Fig. 7. Results of experiments 
are also shown. Each glulam specimen broke at the finger 
joint and the results for other glulam grades were similar. 
The simulation method is that described above. Five hundred 
virtual glulam specimens were generated. The fracture cri¬ 
teria for bending and axial strength are a ni lf ni + oJf bi > 1 and 




Fig. 6. Results of simulation and experiments (type A and B speci¬ 
mens). Solid line, MOR (simulation); circles, MOR (experiment); 
dotted line, <r T (simulation); squares, cr T (experiment); dashed line, a c 
(simulation); crosses, c c (experiment) 
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Fig. 7. Results of simulation and experiments (type C and D 
specimens) 


aJf ni >1, respectively. Here, a m and a,„ are the axial stress 
and the bending stress of the i-th lamina, respectively, and 
f m and f hi are the axial strength and the bending strength of 
the i-th lamina, respectively. Note that the cumulative prob¬ 
ability of the experiments is plotted for convenience, but is 
not particularly meaningful. It can be seen that the experi¬ 
mental measurements are roughly in the same range as the 


simulation, except for compression strengths of type C and 
D specimens. The reason that the simulated values are 
smaller than the experimental values for compression 
strength is not clear, but the fracture criteria for compres¬ 
sion strength might not be appropriate. 


Conclusions 

A simulation method for generating the strength of glulam 
using correlated random variables is described here; the 
correlated random variables are MOE, MOR, <7 T , and cr c , 
among others. Experiments to validate the method were 
conducted. By comparing the simulation and experiment 
results, this method was shown to be useful for generating 
the strength of glulam. 
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